Collective oscillations of a trapped Fermi gas near a Feshbach resonance 
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The frequencies of the coUective oscillations of a harmonically trapped Fermi gas interacting with 
large scattering lengths are calculated at zero temperature using hydrodynamic theory. Different 
regimes are considered, including the molecular Bose-Einstein condensate and the unitarity limit 
for collisions. We show that the frequency of the radial compressional mode in an elongated trap 
exhibits a pronounced non monotonous dependence on the scattering length, reflecting the role of 
the interactions in the equation of state. 

^ ' PACS numbers: 3.75.-b, 3.75.Ss, 67.40 Db 
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The recent realization Q, |^ of a Bose-Einstein condensate of pairs of Fermi atoms (molecular BEC) near a 
^ Fesbach resonance is opening new exciting perspectives in the investigation of highly correlated many-body systems 
^ including the investigation of the long sought BCS-BEC crossover ^4]. An important feature of these configurations is 
the hydrodynamic nature of their macroscopic dynamics which, at zero temperature, is a general consequence of 
superfluidity. Though the observation of the hydrodynamic behaviour cannot be used as a direct proof of superfluidity 
! 0, its careful investigation can provide unique information on the equation of state, thereby checking the predictions 
of many-body theories beyond the mean field picture. 
y—i ' The purpose of this letter is to show that the study of the collective oscillations in a Fermi gas close to a Feshbach 
^ ' resonance is well suited to this purpose, their frequencies exhibiting a non trivial dependence on the scattering length 
\l in regimes now accessible experimentally 0- 

Let us first consider the case of large and positive atomic scattering lengths close to a Feshbach resonance. This 
regime is associated with the occurrence of a weakly bound molecular state and consequently, for sufficently low 
temperatures and low densities, one expects the formation of a Bose-Einstein condensate of molecules. If the molecular 
' gas parameter nMa\i is sufficiently small {ni\i and ai\i are, respectively, the molecular density and the molecule- 
> molecule scattering length) the molecular Bose-Einstein condensate is described by Gross-Pitaevskii theory. In the 
presence of harmonic trapping one finds that in the Thomas-Fermi limit the frequencies of the collective oscillations 
do not depend on the interaction coupling constant, but only on the value of the axial and radial trapping frequencies 
ujz and uj± 0. In the following we will discuss the case of elongated configurations (ujz <SC w±) which are easily 
■ produced experimentally and where one can obtain simple analytic results at zero temperature. The coupled monopole- 
^ ' quadrupole m = modes (m is the third component of angular momentum) oscillate with frequencies 8] corad = 
O • and Uax = '\/5/2wz along the radial and axial directions resp ectively. These modes have been already observed in 
' traditional Bose-Einstein condensates (see, for example, [qLIIOII). First measurements of cOax in a molecular condensate 
^ . have been also reported 0. When the molecular gas parameter tt-mo'^/ is no longer small the above results are modified 
^ by the presence of beyond mean field effects in the equation of state accounted for by the expansion 

........... 32 ^r— 



(N 



fJ-M = gnnM (^1 + ■^-j=\J nMa\i j (1) 

for the molecular chemical potential. Here — 2}T?aM/M is the molecule-molecule interaction coupling constant 
and M is the atomic mass. In terms of the atomic chemical potential /i and density n one has /ijv/ — 2fi and um — n/2. 
The relation between the molecule-molecule scattering length aj\/ and the atomic scattering length a has been recently 
calculated in j3| smd is given by um — 0.6a. The collective frequencies predicted by the new equation of state in 
the presence of harmonic trapping can be calculated using hydrodynamic theory. The result for the m — radial 
frequency is 
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The relative correction to the mean field value \fhj2LOz of the to = axial frequency u^ax turns out to be a factor 6 
smaller, suggesting that the radial mode is more suited to explore the interaction effects in the equation of state [l^ . 
In the above equation nM(0) is the density of the molecular condensate evaluated in the center of the trap. It is useful 
to express the molecular gas parameter in terms of the relevant atomic parameters. By using the Thomas-Fermi value 
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for the central density one finds 
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nM(0)ai, ^O.Oe hvi/^— (3) 

where N is the total number of atoms and aho = \/hjMUj is the atomic oscillator length relative to the geometrical 
average w of the trapping frequencies. When the scattering length a is further increased and becomes comparable or 
even larger than the average distance d ^ n~^/^ between particles the molecular description is no longer applicable. 
Eventually, when a ^ d the system reaches the so called unitarity limit for collisions where the equation of state does 
not depend any more on the actual value of the scattering length, nor on its sign. Provided the other relevant lengths 
associated with the two-body interaction, like the range of the potential, are smaller than d, the equation of state is 
expected to exhibit a universal behaviour Q . Considerations of dimensionality suggest that at zero temperature the 
chemical potential should have the same density dependence 

A^=(l + /3)|^(6^'f/'n2/3 (4) 

as in the ideal Fermi gas, the parameter /? being a universal coefficient accounting for quantum and dynamical 
correlations. At present the most reliable theoretical estimate of (3 for a two component Fermi system is /? ^ —0.56 
obtained in through a correlated basis function calculation. Actually the value of /3 is not relevant for the 
determination of the collective frequency 'l6] which depends only on the power law coefficient characterizing the 
density dependence of the chemical potential. Ifu oc ii?/^ hydrodynamic theory predicts, for an elongated trap, the 
values L Orad = \/ 10 /3uj± and ujax = v^lS/Sw^ [131 . The same results hold for a classical gas in the hydrodynamic 
regime where, at constant entropy, the chemical potential is also proportional to n^/^. 

It is worth noticing that the value predicted for ujrad is lower than the value 2uj± holding in the molecular BEC 
regime, revealing that the frequency of the radial compression mode first increases as a function of a, as predicted by 
(0) , and then decreases to reach the value y/l0/3u!± in the unitarity limit (see the Figure) . 

Let us finally discuss the region of negative scattering lengths, corresponding to the BCS regime of superfluidity. 
When the gas parameter n \ a ^ is much smaller than unity the chemical potential still exhibits a 2/3 power law 
dependence on the density and one consequently expects to find again the result uJrad — \/W/3uj± in the hydrodynamic 
regime. However, differently from the case of positive scattering lengths, in this limit the critical temperature for 
superfluidity becomes very small, and the superfluid phase is likely much more difficult to reach. Mean field corrections 
to the equation of state can be also easily calculated and yield the expansion 

^^=^i6nY/'ny^ + \9n (5) 

for the chemical potential where g = ATrfi^a/M. By treating the interaction as a small perturbation one finds the 
following result for the frequency of the radial compression mode: 

where Eint — (ff/4) j dru? and Eho — J drVhon are, respectively, the interaction and oscillator energies. In terms of 
the dimensionless combination N^^^a/uho one has Eint/ E^o = 0.5N^^^a/aho The mean field correction in lO is 
a factor 4/5 smaller than the analogous term calculated for the monopole mode in the case of isotropic trapping (20j. 
The interaction correction in curad is negative, reflecting the attractive nature of the force, so that ujrad first decreases 
as a function of | a | and then increases to reach the value y^l0/3u!± in the unitarity limit (see the Figure) . 

The determination of tOrad in the intermediate regimes requires the full knowledge of the equation of state, through 
a consistent many-body calculation. In the Figure we have drawn a schematic interpolation (dot-dashed line) between 
the asymptotic behaviours Q and © through the resonance. The curve explicitly points out the oscillatory behaviour 
of the frequency. 

In conclusion we have shown that the frequency of the radial compression mode in an ultracold trapped Fermi gas 
exhibits a non trivial dependence on the scattering length near a Feshbach resonance, reflecting the occurrence of 
important interaction effects in the equation of state. These effects characterize regimes of high physical interest, now 
available experimentally. 

Stimulating discussions with R. Grimm are acknowledged. 
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FIG. 1: Radial compressional frequency in an elongated trapped Fermi gas close to a Feshbach resonance as a function of the 
dimensionless parameter {N^^^a/aho)~^ ■ The mean field theory in the molecular BEG regime predicts uj/lu± = 2, while in the 
unitarity limit one expects lu/uj± — \/10/3 ~ 1.83. The full line for a > corresponds to the expansion Q accounting for 
beyond mean field (m.f.) corrections in the molecular condensate with sa/ ~ 0.6a, while the full line for a < corresponds to 
the expansion Q accounting for the mean field correction in the BGS phase. The dot-dashed line is a schematic interpolation 
between the two asymptotic regimes. 
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